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1.1 Introduction 

The scientific belief that the universe evolves in time is one of the legacies 
of the theory of the Big Bang. 

The concept that the universe has an history started to attract the 
interest of cosmologists soon after the first formulation of the theory: al- 
ready Gamow (1948; 1949) investigated how and when galaxies could have 
been formed in the context of the expanding Universe. 

However, the specific topic of the formation (and of the fate) of the first 
objects dates to two decades later, when no objects with metallicities as low 
as those predicted by primordial nucleosynthesis {Z ^ 10~^° ~ 10~^Zq) 
were found. Such concerns were addressed in two seminal papers by Pee- 
bles & Dicke (1968; hereafter PD68) and by Doroshkevich, Zel'Dovich & 
Novikov (1967; hereafter DZN67)^, introducing the idea that some objects 
could have formed before the stars we presently observe. 

^ This paper is in Russian and we base our comments on indirect Icnowledge {e.g. from 
the Novikov & Zel'Dovich 1967 review). 
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(i) Both PD68 and DZN67 suggest a mass of ~ lO^M© for the first genera- 
tion of bound systems, based on the considerations on the cosmological 
Jeans length (Gamow 1948; Peebles 1965) and the possible shape of 
the power spectrum. 

(ii) They point out the role of thermal instabilities in the formation of 
the proto-galactic bound object, and of the cooling of the gas inside 
it; in particular, PD68 introduces H2 cooling and chemistry in the 
calculations about the contraction of the gas. 

(iii) Even if they do not specifically address the occurrence of fragmenta- 
tion, these papers make two very different assumptions: PD68 assumes 
that the gas will fragment into "normal" stars to form globular clus- 
ters, while DZN67 assumes that fragmentation does not occur, and 
that a single "super-star" forms. 

(iv) Finally, some feedback effects as considered {e.g. Peebles & Dicke con- 
sidered the effects of supernovae). 



Today most of the research focuses on the issues when fragmentation 
may occur, what objects are formed and how they influence subsequent 
structure formation. 

In these notes we will leave the discussion of feedback to lecture notes 

by Ferrara & Salvaterra and by Madau & Haardt in this same book and 
focus only on the aspects of the formation of the first objects. 

The advent of cosmological numerical hydrodynamics in particular al- 
low a fresh new look at these questions. Hence, these notes will touch on 
aspects of theoretical cosmology to chemistry, computer science, hydrody- 
namics and atomic physics. For further reading and more references on 
the subject we refer the reader to other relevant reviews such as Barkana 
& Loeb 2001, and more recently Ciardi & Ferrara 2004, Glover 2004 and 
Bromm & Larson 2004. 

In this notes, we try to give a brief introduction to only the most 
relevant aspects. We will start with a brief overview of the relevant cosmo- 
logical concepts in section 2, followed by a discussion of the properties of 
primordial material (with particular emphasis to its cooling and its chem- 
istry) in section 3. We will then review the technique and the results of 
numerical simulations in sections 4 and 5: the former will deal with detailed 
3D simulations of the formation of gaseous clouds which are likely to trans- 
form into luminous objects, while the latter will examine results (mostly 
from ID codes) about the modalities of such transformation. Finally, in 
section 6 we will critically discuss the results of the previous sections, ex- 
amining their consequences and comparing them to our present knowledge 
of the universe. 
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1.2 Physical cosmology 

In the following of this notes we will adopt the modern physical description 
of the universe (dating back to at least 1920), based upon the "cosmological 
principle" , which affirms that on cosmological scales the distributions of 
matter and energy should be homogeneous and isotropic, whose metric is 
the Robertson- Walker metric. Although derived from mainly philosophical 
arguments, the cosmological principle is also supported by observations 
such as the isotropy of the CMB {e.g. Wu et al. 1999). 

We will also make some additional, general assumptions which are 
quite common in present-day cosmology, and which are believed to be the 
best explanations for a scries of observations. That is: 

(i) The cosmological structures we observe at present (galaxies, clusters 
of galaxies etc.) formed because of gravitational instability of pre- 
existing, much shallower fluctuations; 

(ii) Most of the matter in the universe is in the form of "Cold Dark Mat- 
ter", that is, of some kind of elementary particle (or particles) that 
has not been discovered at present. Cold Dark Matter particles are 
assumed to have a negligibly small cross section for electromagnetic 
interactions {i.e. to be dark), and to move at non-relativistic speeds 
{i.e. to hecold). 

1.2.1 Fluctuations in the Early Universe 

1.2.1.1 Inflation 

Inflation is a mechanism which was first proposed by Guth (1981) and (in 
a different way) by Starobinsky (1979, 1980), and has since been included 
in a number of different "inflationary theories" (see Guth 2004 for a review 
upon which we base this paragraph). 

The basic assumption of inflationary models is the existence of states 
with negative pressure; a typical explanation is that some unidentified kind 
of scalar field (commonly referred to as inflaton) temporarily keeps part of 
the universe in a "false vacuum" state, in which the energy density must 
be approximately constant (that is, it can not "diluted" by expansion) at 
some value p/c^, implying a negative pressure pf = —pfC^. Inserting these 
two expression in the first Friedmann cosmological equation 

a{t) = -fG(p + 3^)a{t) (1.1) 

(where a{t) is the scale factor at cosmic time; t) it is easy to derive that the 
considered region expands exponentially: a{t) oc e*/'-' withtf = {SnGpf /3)~^^^; 
the epoch in which the universe undergoes such exponential expansion is 
called inflation. 



Physical cosmology 4 




Figure 1.1. Effect of universe geometry on tlie observed angular size of fluc- 
tuations in tlie CMBR. If tiie universe is closed (left panel) "hot spots" appear 
larger than actual size; if the universe is flat (middle panel), "hot spots" appear 
with their actual size; if the universe is open, "hot spots" appear smaller than 
their actual size. 

If p/ is taken to be at about the grand unified theory scale, we get 
tf ~ 10""^* s, corresponding to an Hubble length of ct/ ~ 10~^^ cm; 
if the inflationary phase is long enough (a lower limit for this is about 
65 e-folding times, corresponding to an expansion by a factor ~ 10^*), it 
smooths the metric, so that the expanding region approaches a de Sitter flat 
space, regardless of initial conditions. Furthermore, when inflation ends, 
the energy stored in the inflaton field is finally released, thermalizes and 
leads to the hot mixture of particles assumed by the standard big bang 
picture. 

Inflation helps explaining several fundamental observations, which were 
just assumed as initial conditions in non inflationary models: 

(i) The Hubble expansion: repulsive gravity associated with false vacuum 
is exactly the kind of force needed to set up a motion pattern in which 
every two particles are moving apart with a velocity proportional to 
their distance. 

(ii) The homogeneity and isotropy: in "classical" cosmology the remark- 
able uniformity of the Cosmic Microwave Background Radiation (CMBR) 
cannot be "explained" , because different regions of the present CMBR 
sky never were causally connected. Instead, in inflationary models the 
whole CMBR sky was causally connected before inflation started, and 
uniformity can be established at that time. 
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(iii) The flatness problem: a flat Friedman-Robertson- Walker model uni- 
verse (i.e. with Cl(t) = ptot{t) / Pcritit) — I7 where ptot{t) is the cosmo- 
logical mean density, including the "dark energy" term, and Pcrit(^) = 
3H{t)'^ /SttG, H{t) = a{t)/a{t) being the Hubble parameter) always 
remains flat, but if at an early time 17 was just slightly different from 
1, the difference should have grown as fast as (17 - 1) oc t^/s. aU 
the observational results {e.g. the Bennet et al. 2003 WMAP result 
i7o = 1.02 ± 0.02) show that at present is quite close to 1, imply- 
ing that at the Planck time {t ~ 10"^'^ s) the closeness was really 
amazing. Inflation removes the need for this incredibly acctirate flne 
tuning, since during an inflationary phase Q. is driven towards 1 as 
(17 - 1) oc e-^"K 

(iv) The absence of magnetic mono-poles: grand unified theories, combined 

with classical cosmology, predict that the universe should be largely 
dominated by magnetic mono-poles, which instead have never been 
observed. Inflation provides an explanation, since it can dilute the 
mono-pole density to completely negligible levels. 

(v) the anisotropy properties of the CMBR radiation: inflation provides a 
prediction for the power spectrum of fluctuations, which should be gen- 
erated by quantum fluctuations and nearly scale-invariant, a prediction 
in close agreement with the WMAP results (see the next subsection). 

A peculiar property of inflation is that most inflationary models pre- 
dict that inflation does not stop everywhere at the same time, but just in 
localized "patches" in a succession which continues eternally; since each 
"patch" (such as the one we would be living in) can be considered a whole 
universe, it can be said that inflation produces an infinite number of uni- 
verses. 



1.2.1.2 Primordial fluctuation evolution - Dark Matter 

Inflation predicts the statistical properties of the density perturbation field, 
defined as 

j(x) ^ eiil-I (1.2) 

where r is the; proper coordinate, x = r/a is the comoving coordinate, and 
p is the mean matter density. 

In fact, if we look at the Fourier components in terms of the comoving 
wave-vectors k 

^(k) = J 5(x) e"'""" rf^x (1.3) 

the inflationary prediction is that the perturbation field is a Gaussian ran- 
dom field, that the various k modes are independent, and that the power 
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Figure 1.2. Schematic siiape of tlie DM power spectrum after accounting for 
tiie effects of free streaming and of the processing that takes place for fluctuations 
entering the horizon before icq- 



spectrum P{k) (where k = |k|) is close to scale-invariant, i.e., it is given 
by a power law 

P{k) =< |(5kP >oc fc"% with ris ~ 1. (1.4) 

This prediction applies to the primordial power spectrum; in order 
to make comparisons with observations, we need to include effects which 
subsequently modified its shape: 

(i) Free streaming: the dark matter particles are in motion, and since 
they are believed to interact only through gravity, they freely propa- 
gate from overdense to underdense regions, wiping out the density per- 
turbations. However, at any given cosmic time t, this will affect only 
wavelengths smaller than the free streaming length, i.e. the proper 
distance ^ps (t) which a dark matter particle can have travelled in time 
t. It can be shown (see e.g. Padmanabhan 1993, §4.6) that this scale 
depends primarily on the epoch t^r when the dark matter particles 
become non-relativistic: before that epoch they move at the speed of 
light, and can cover a proper distance ZFs(^nr) ~ Scim-, correspond- 
ing to a comoving distance Aps = (ao/flnr)2cinr; after becoming non- 
relativistic, particle motions become much slower than cosmological 
expansion, and the evolution after tnr only increases Aps by a rela- 
tively small factor. In turn, t^^ depends primarily on the mass todm 
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of dark matter particles, so that 

~ 5 X 10-3 {n^uhy^^ (r^)"'^' (1-5) 

corresponding to a mass scale 

Mfs 6 X 10-15 (Qdm/i') (f^) ' Mq. (1.6) 

Since the most favoured candidates for the DM particles arc Weakly 
Interacting Massive Particles (WIMPs) with mass between 0.1 and 
100 GeV, we probably have Mfs rS 10~^ Mq. Some super-symmetric 
theories {e.g. Schwartz et al. 2001, and more recently Green, Hof- 
mann & Schwarz 2004) instead point towards Mps ~ 1Q~^Mq, and 
Diemand, Moore & Stadel (2005) used this result in order to argue that 
the first structure which formed in the early Universe were Earth-mass 
dark-matter haloes (but see also Zhao et al. 2005 for criticism about 
this result); 

(ii) Growth rate changes: in general, perturbations grow because of grav- 
ity; however, the details of the process change with time, leaving an 
imprint on the final processed spectrum. The time teq when the mat- 
ter density becomes larger than the radiation density is particularly 
important: before feqi perturbations with A larger than the Hub- 
ble radius c/H{t) grow as ^ oc a^, while the growth of smaller per- 
turbations is almost completely suppressed; after tcq both large and 
small fluctuations grow as S (x a. Because of these differences, the 
size of the Hubble radius at tgq, Xeq — 13 (f2/i^)-^Mpc (in terms of 
mass, Mcq ~ 3.2 x 10"'^^(f7/i^)-^ Mq) separates two different spectral 
regimes. In the wavelength range Afs < A < Aeq the growth of fluc- 
tuations pauses between the time they enter the Hubble radius and 
teq. As a result the slope of the processed spectrum is changed, and 
P{k) oc A;"=— * ~ k~^. Instead, at scales A > Aeq all the fluctuations 
keep growing at the same rate at all times, and the shape of the power 
spectrum remains similar to the primordial one, P{k) oc k"" ~ k^. 

WMAP (Bennett et al. 2003) measured the spectral index, obtaining 
Hs = 0.99 ± 0.04, and did not detect deviations from gaussianity, both 
results in agreement with inflationary predictions. 

This kind of spectrum, in which fluctuations are typic:ally larger on 
small scales, leads naturally to hierarchical structure formation, since small- 
scale fluctuations are the first to become non-linear {i.e. to reach 6 ~ 
1), collapse and form some kind of astronomical object. It is also worth 
remarking that the very first objects, coming from the highest peaks of 
(5(x), are typically located where modes i5(k) of different wavelength make 
some kind of "constructive interference" : the very first objects are likely to 
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Figure 1.3. Map of the temperature anisotropies in tlie CMB as obtained by 
combining the five bands of the WMAP satellite in order to minimise foreground 
contamination (from Bennett et al. 2003). 

be on top of larger ones, and they are likely to be clustered together, rather 
than uniformly distributed. For this reason, it is also very likely that the 
halos where these objects formed have long since been incorporated inside 
larger objects, such as the bulges of Af, galaxies or the cD galaxy at the 
centre of galaxy clusters (see e.g. White & Springel 1999). 



1.2.1.3 Fluctuation evolution - Baryons 

Before the equivalence epoch icq the baryons evolve in the same way as 
dark matter. Instead, in the matter dominated era they behave differently: 
we mentioned that all the dark matter fluctuations which were not erased 
by free streaming grow as 6 (x a oc [1 + z)~^ , but this does not apply to 
baryons. In fact, baryons decouple from radiation only at <doc, significantly 
later than tcq (we remind that 1 + Zcq ^ 10^, while 1 + z^cc — lO'^). The 
persistence of the coupling with radiation prevents the growth of baryonic 
fluctuations on all scales; even worse, on relatively small scales all the 
fluctuations in the baryonic component are erased through a mechanism 
similar to free streaming. Such effect takes place below the so-called Silk 
scale (Silk 1968), which is given by the average distance that the photons 
(and the baryons coupled with them) can diffuse before t = idee! this 
translates into a comoving distance 



As ^ 3.5 




(1.7) 
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Figure 1.4. Angular CMB power spectrum of temperature (top panel) and 
temperature-polarization (bottom panel) as obtained by the WMAP satellite 
(from Bennett et al. 2003). The line shows the best-fit with a ACDM model, and 
grey areas represent cosmic variance. 



(where fib is the baryonic contribution to fi) and encloses a mass 



Ms ~ 6.2 X 10^^ 



n 



3/2 



(1.8) 



This result was a major problem for cosmology before the existence 
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of dark matter started to be assumed: it implies that in a purely baryonic 
universe there should be no structures on a scale smaller than that of 
galaxy clusters (if f2 = fib ~ 0.1). Furthermore, even fluctuations which 
were not erased can grow only by a factor 1 + 2;dec between decoupling and 
present time, and this is not enough to take a typical CMBR fluctuation (of 
amplitude 5 ~ 10~^) into the non-linear regime. The introduction of Cold 
Dark matter solved this problem, since after recombination the baryons 
are finally free to fall inside dark matter potential wells, whose growth was 
unaffected by the photon drag. 

It can be found that after decoupling from radiation, the baryonic 
fluctuations quickly "reach" the levels of dark matter fluctuations, evolving 
as 

6b = Sum f 1 - t^^) , (1-9) 

so that the existing dark matter potential wells "regenerate" baryonic fluc- 
tuations, including the ones below the Silk scale. 

This result is correct as long as pressure does not play a role, that 
is for objects with a mass larger than the cosmological Jeans mass Mj oc 
7^3/2^-1/2 _ g^j(.jj mass behaves differently at high and low redshift. Before 
a redshift ^Jcompton — 130 we have that the temperature of the baryons is 
still coupled to that of the CMB because of Compton scattering of radiation 
on the residual free electrons; for this reason, Tb{z) ~ Tcmb{z) oc (1 + z), 
and as p{z) oc (1 -|- z)^ the value of Mj is constant: 

Mj(^) ~ 1.35 X 10^ (^^^^ ^ Mq (for ^dec ^ ^ ^ ^Compton) 

(1.10) 

where fim = fib + ^^dm is the total matter density (baryons plus dark 
matter). At lower redshifts the baryon temperature is no longer locked to 
that of the CMB and drops adiabatically as Tb oc (1-1-2:)^. At such redshifts 
the Jeans mass evolves as 

„,,.,..„.(l^)- (-!)-(-)'%,,,„, 

(for Z < ZCompton). 



1.2.2 From fluctuations to cosmological structures 

1.2.2.1 Non-linear evolution 

When the density contrast S becomes comparable to unity, the evolution 
becomes non-linear, and the Fourier modes no more evolve independently. 
The easiest way to study this phase is through the "real space" d{x) (rather 
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Figure 1.5. Schematic evolution of the temperature of the Intergalactic Medium 
(IGM) as a function of redshift: after recombination and before thermal decou- 
pling (at zcompton), TiGM is locked to Tcmbr and evolves as (l + z). After thermal 
decoupling Tigm is no more locked to the radiation temperature, and drops adi- 
abatically as (1 + z)^, until the first non-linear objects collapse, and emit light 
which is able to re-ionize the universe, raising Tigm above ~ 10* K. 



than its Fourier components 5(k)), considering the idealized case of the col- 
lapse of a sphericaUy symmetric overdensity, and in particular the collapse 
of a top-hat fluctuation. It is well known that, through some simple further 
assumption (such as that "spherical shells do not cross" ) , the time evolu- 
tion of the radius i? of a top-hat perturbation (see e.g. Padmabhan 1993, 
§8.2 for the treatment of a slightly more general case) can be written down 
as 

= !o-i' ,A ^-^o,9it)] (1.12) 

2 6^- {a, - 1) 

where Ri = R{ti), Si = 5{ti) and Jl^ = ri(ti) are the "initial conditions" for 
the fluctuation evolution, and 9 is defined by the equation 



[e{t) - sinOit)]- 



1 + 5, 



t 



(1.13) 



where again Hi is the Hubble parameter at ti, and the last approximate 
equality is valid only as long as 5i <C 1 (that is, a sufficiently early ti must 
be chosen). The fluctuation radius R reaches a maximum _Rta at the so- 
called turn-around epoch (when 9 = n) when the overdense region finally 
detaches itself from the Hubble flow and starts contracting. However, while 
eq. (|1.12|l suggests an infinite contraction to i? = (when 9 = 2tt), the vi- 
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olent relaxation process (Lynden-Bell 1967) prevents this from happening, 
leading to a configuration in virial equilibrium at i?vir — 

Here, we summarise some well known, useful findings of this model. 

First of all, combining the evolution of the background density evo- 
lution and of eq. (|1.12() it is possible to estimate the density contrast 
evolution 

which leads to some noteworthy observation, such as that the density con- 
trast at turn-around is Sta. = (97r^/16) — 1 ~ 4.6, which at virialization 
becomes S^ir = Ac, where it can be usually assumed that Ac ~ ISvr^ but 
sometimes higher order approximations are necessary, such as the one in 
Bryan & Norman 1998, 

Ac = IStt^ + 82(1 - n^^) ~ 39(1 - n'J^ (1.15) 

with 

where i^A is the dark energy density, and flk — 1 — — J7a is the curvature. 
From this, it is possible to estimate the virial radius 

(1.17) 

the circular velocity for such an halo 

and the virial temperature 
T., . 19800 { f^^Y' fi±i) f^) K. (1.19) 

1.2.2.2 The Press- Schechter formalism 

The simple top-hat model is at the core of the so-called Press-Schechter 
formalism (Press & Schechter 1974, but see also the contribution by Som- 
merville in this same book), which predicts the density of virialized halos 
of a given mass at a given redshift. This model assumes that the distribu- 
tion of the smoothed density field Sm (where M is the mass scale of the 
smoothing) at a certain redshift zq is Gaussian with a variance ctai , so that 
the probability of having S m larger than a given Smt is 



(1.18) 
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Figure 1.6. Characteristic mass of la (bottom solid curve), 2a (middle solid 
curve) and 3a (top solid curve) collapsing halos as a function of redshift. These 
were obtained from the Eisenstein & Hu (1999) power spectrum, assuming an 
Qa = 0.7, rim = 0.3 cosmology. The dashed curves show the minimum mass 
which is required for the baryons to be able to cool and collapse (see next section) 
in case of pure atomic cooling (upper curve) and of molecular cooling (lower 
curve), (from Barkana & Loeb 2001). 



a common choice is zq — 0, requiring Sm to be estimated through a purely 
linear evolution of the primordial power spectrum. 

The Press-Schechter model then chooses a Jcrit = ^crit(-z) (but it is also 
possible to assume a constant Sent and make cta/ a function of redshift; see 
e.g. Viana & Liddle 1996) and assumes that this probability (multiplied by 
a factor of 2 - see Bond et al. 1991 for an explanation of this extra factor) 
also gives the fraction mass which at a redshift z is inside virialized halos 
of mass M or larger. This can be differentiated over M in order to get the 
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Figure 1.7. Halo mass functions at several redshifts (from bottom left to bottom 
right, z = 30, z — 20, z — 10, z — 5 and z = 0, respectively). The assumed power 
spectrum and cosmology are the same as in 

fig- EH (from Barkana & Loeb 2001). 

mass distribution at each given redshift 

dn _ 2 p„i dln{l/aM) ScTitjz) - ^""a^''' n 211 

dM ~ (27r)i/2 M dM am ^ " 

In this way, the abundance of halos is completely determined through 
the two functions Scrit{z) and ctm- The first one is commonly written as 
i5crit(.2) — So/D{z), where D{z) is the growth factor {D{z) ~ (1 + z)^^ for 
Einstein-de Sitter models; see Peebles 1993 for a more general expression), 
coming from cosmology; instead is usually taken to be 1.686, since the 
top-hat model predicts that an object virializes at the time when the linear 
theory estimates its overdensity at S = 1.686. Instead ctm depends from 
the power spectrum; for example, figures 11.61 and 11.71 are based on the 
Eisenstein & Hu (1999) results^. 

^ The authors of this paper also provide some very useful codes for dealing 
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1.3 Primordial gas properties 

1.3.1 Cooling 

The typical densities reached by the gas after viriaUzation (of the order of 
riB = ps/mii ~ O.Oirib[(l + 2;vir)/10]^cm~^ are far too low for the gas to 
condense and form an object like a star. The only way to proceed further 
in the collapse and in the formation of luminous objects is to remove the 
gas thermal energy through radiative cooling. 

For this reason, cooling processes are important in determining where, 
when and how the first luminous objects will form. 

In Fig. II. Sl it is possible to see that the cooling of primordial {i.e. metal- 
free) atomic gas at temperatures below 10^ K is dramatically inefficient, 
because in that temperature range the gap between the fundamental and 
the lowest excited levels (~ 10.2 eV for H atoms) is so much larger than 
the thermal energy ~ A:bT ^ 1 eV that very few atoms get collisionally 
excited, and very few photons are emitted through the corresponding de- 
excitations. 

This is important because in all hierarchical scenarios the first objects 
to virialize are the smallest ones, and such halos have the lowest virial tem- 
peratures. If the primordial gas were completely atomic, the first luminous 
objects would probably form relatively late (z ^ 10 — 15), in moderately 
massive halos (Af ~ 10^ Mq) with T^i, > 10^ K. 

However, it is also possible to see from Fig. 11.81 that the presence 
of molecules in small amounts (/hj = 2nH2/("H + SnRa) ^ 5 x 10~^; the 
dashed curve in Fig. II. 81 was obtained assuming fn^ — 10^^) can dramat- 
ically affect the cooling properties of primordial gas at low temperatures, 
making low mass halos virializing at high redshift (z ^ 20) the most likely 
sites for the formation of the first luminous objects. 

1.3.2 Molecular cooling 

In the current scenario for the formation of primordial objects, the most 
relevant molecule is H2. The only reason for this is the high abundance of 
H2 when compared with all other molecules. In fact, the radiating prop- 
erties of an II2 molecule are very poor: because of the absence of a dipole 
moment, radiation is emitted only through weak quadrupole transitions. In 
addition, the energy difference between the H2 ground state and the lowest 
H2 excited roto- vibrational levels is relatively large (AiJoi/fce ^ 200 K, be- 
tween the fundamental and the first excited level; however, such transition 
is prohibited by quadrupole selection rules, and the lowest energy gap for 
a quadrupole transition is A£^o2/^b — 510 K), further reducing the cooling 
efficiency at low temperatures. 

with the power spectrum and the Press-Schechter formalism at the web page 
|http://background. uchicago.edu/~whu/transfer/transferpage. html] 
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Figure 1.8. Cooling rate per atomic mass unit of metal-free gas, as a function of 
temperature. The solid line shows assumes the gas to be completely atomic (the 
two peaks correspond to H and He excitations, while the high temperature tail is 
dominated by free-free processes) and drops to about zero below T ~ 10'' K; the 
dashed line shows the contribution of a small (/h2 = 10""^) fraction of molecular 
hydrogen, which contributes extra-cooling in the range 100 K ^ T ^ 10^ K (from 
Barkana & Loeb 2001). 



Apart from H2 the most relevant molecular coolants are HD and LiH. 
In the following, we will briefly list the cooling rates (mainly taken from 
GalH & Palla 1998, hereafter GP98; see also HoUenbach & McKee 1979, 
Lepp & ShuU 1983, 1984 and Martin et al. 1996, Le Bourlot et al. 1999, 
Flower et al. 2000) of H2 and of the other two possibly relevant species'^. 

We note that Bromm et al. 1999 included HD cooling in some of their 

^ In Fig. II.IUI the cooling rate for H^is shown, too. But while it is still marginally 
possible that HD or LiH cooling can play some kind of role in the primordial universe, 
this is much more unlikely for H J , mainly because its under-abundance with respect to 
H2 is always much larger than the difference in the cooling rates; for this reasons we 
choose to omit a detailed discussion of H J cooling. 
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Lepp and Shull (1983, red), Martin et al. (1996, blue), Galli and Palla (1998, green) 
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Figure 1.9. Cooling rate per H2 molecule as computed by different authors 
(Lepp & Shull 1983, Martin et al. 1996, GP98) as a function of temperature 



and for different densities; note that for n ^ 
independent of density. 
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simulations but found that it never accounted for more than ~ 10% of 
H2 cooling; however, they did not completely rule out the possibility that 
a "HD /'^H2 I'atio substantially larger than the equilibrium value (close to 
nD/nn ~ 10~^) could change this conclusion. Also, in present-day simula- 
tions of primordial star formation the gas temperature never goes below a 
few hundreds degrees: in case it did, even a tiny amount of LiH could be 
enough to dominate the cooling rate. 

1.3.2.1 H2 cooling rate 

The H2 cooling rate per molecule (p, T) can be conveniently expressed 
in the form: 

AH.(P,7^) = -^^fg^ (1.22) 

where Ah2,lte(T') and nuAHa.p-^o are the high and low density limits of 
the cooling rate (which apply at n 10'' cm^'^ and at n ;$ 10^ cm~^, re- 
spectively). 
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The high density (or LTE, from the Local Thermal Equilibrium as- 
sumption which holds in these conditions) limit of the cooling rate per 
H2 molecule is given by HoUenbach & McKee (1979): 

OA 

^ ergs"^ (1.23) 

where T3 = T/(1000 K). Note that the first row in the formula accounts 
for rotational cooling, while the second row accounts for the first two vi- 
brational terms. 

For the low density limit, GP98 found that in the relevant temperature 
range (10 K < T < 10^ K) the cooling rate Ana.p^o is independent from 
density, and is well approximated by 

logAH2,p^o(T) ~ -103.0 97.59 logT - 48.05(logT)2 -h 

10.80(logT)^ - 0.9032(logT)^ (1.24) 

where T and Ah2,p-»o are expressed in K and ergs^^ cm'', respectively. 

Note that even if both AhjXTE and Ahs.p^o do not depend on density, 
Ah2(p, T) is independent of p only in the high density limit. 

1.3.2.2 HD and LiH cooling rates 

The cooling rates of HD and of LiH are more complicated (see Flower et 
al. 2000 for HD and Bogleux & GaUi 1997 for LiH), but in the low density 
limit (and in the temperature range 10 K < T < 1000 K) it is possible to 
use the relatively simple expressions given by GP98. 

For HD, we have that the low density limit of the cooling rate per 
molecule, Ahd.p^O! is: 

Ahd,p^o(T) ~ 27io£;ioe"^ + {5/3)j2iE2ie~^ (1.25) 

where E'lo and E21 are the energy gaps between HD levels 1 and and 
levels 2 and 1, respectively; they are usually expressed as Eio = ksTio and 
E21 = A:bT2i, with Tio 128 K and T21 ~ 255 K, 

710 and 721 are the approximate coUisional de-excitation rates for the 
1-0 and 2-1 transitions, and are given by 

710 4:.4: X 10"^^ + 3.6 X io-13tO-77 (1.26) 
721 ~ 4.1 X 10"^^ + 2.1 X io-13tO-92 (1.27) 

where we use the numerical value of T (in Kelvin) and the rates are ex- 
pressed in cm^ s~^. 



+6.7 X IQ-^^e"^ + 1.6 x lO'^^e^ 
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Figure 1.10. Comparison of the low density (n ^ 10 cm~^) cooling rates per 
molecule of several molecular species, in particular H2 , HD and LiH (from GP98). 
Note that at T 100 K both HD and LiH molecules are more than 10^ times 
more efEcient coolants than H2 molecules, but this difference is believed to be 
compensated by the much higher H2 abundance (see e.g., Bromm et al. 2002). 
The plot also shows the cooling due to H-Hjand e~-Hjcollisions, but these con- 
tributions are never important because of the very low H^abundance. 



For LiH instead we have that the same density limit of the cooling rate 
per molecule, ALiH,p^o can be fitted by: 

logio(ALiH,^,^o) = Co + Ci logio T + C2(logio Tf + 

C3(logioT)3 + C4(logioT)* (1.28) 

where co = -31.47, ci = 8.817, C2 = -4.144, C3 = 0.8292 and a = 
—0.04996, assuming that T is expressed in Kelvin, and that Alih.p^o is 
expressed in ergcm^ s~^. 
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Figure 1.11. CIE cooling for different collisions and different molecular frac- 
tions(0.999, 0.5, 0.001). Top panels show the cooling rates per unit mass: the 
thick solid line is the total CIE cooling, while thin lines show the various compo- 
nents: H2 -If2 (solid), H2 -He (short dashed), H2 -H (dot dashed) and H-He (long 
dashed); the dotted line shows the results of the approximate formula given in 
the text. In the bottom panels the ratios of the various component to H2 -II2 CIE 
is shown. All quantities were calculated assuming n = 10^* cm~^ and X = 0.75 
(from Ripamonti & Abel 2004). 



1.3.2.3 Cooling at high densities: Collision Induced Emission 

During the formation of a protostar an important role is played by the 
so called Collision-Induced Emission (CIE; very often known as Collision- 
Induced Absorption, or CIA), a process which requires pretty high densi- 
ties {n — p/m-n ^ 10^^ — 10^'* cm"'^) to become important (see Lenzuni 
et at 1991, Frommold 1993, Ripamonti & Abel 2004). In fact, Collision- 
Induced Emission takes place when a collision between two H2 molecules 
(or II2 and H, or H2 and He, or even H and He) takes place: during 
the collision, the interacting pair briefly acts as a "super-molecule" with a 
nonzero electric dipole, and a much higher probability of emitting (CIE) or 
absorbing (CIA) a photon than an unperturbed H2 molecule (whose dipole 
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is 0). Because of the very short durations of the colhsions, ( ^ 10~^^ s), 
such mechanism can compete with "normal" emission only in high den- 
sity environments. Furthermore, because of the short durations of the 
interactions, collision-induced lines become very broad and merge into a 
continuum: the H2 CIE spectrum only shows broad peaks corresponding 
to vibrational bands, rather than a large number of narrow roto- vibrational 
lines. This is also important because self-absorption is much less rel- 
evant than for line emission, and in primordial proto-stars simulations 
CIE cooling can be treated with the optically thin approximation up to 
n ~ 10^^ cm^'^ (for H2 lines, the optically thin approximation breaks at 
about n ~ lO^*' cm^^). In figure [TTTl we show the CIE cooling rate for 
a gas with n = 10^^ cm~'^, as a function of temperature and for differ- 
ent chemical compositions; at temperatures between 400 and 7000 K. For 
H2 abundances fn^ = 2nH2/("H+ +"-H + 2nH2) ^0.5 the total CIE cooling 
rate can be approximated by the simple expression (see Ripamonti & Abel 
2004) 

Lcie(p,T,X,/h,) ^ 0.072pT4x/H, ergg-^s"! (1.29) 

where the density is in gcm~'^, the temperature is in K, and X 2± 0.75 is 
the hydrogen fraction (by mass). 

1.3.3 Chemistry 

Since molecules play such an important role in the formation of the first 
luminous objects, it is important to include a proper treatment of their 
abundance evolution. However, a full treatment should keep into account 
~ 20 different species and ^ 100 different reactions. For instance, GP98 
give a chemical network of 87 reactions, which includes e~, H, 11+ , H~, 
D, D+, He, He+, He++, Li, Li+, Li", H2 , H+, HD, HD+, HeH+, LiH, 
LiH+, H+ and H2D+; even their minimal model is too complicated to be 
described here, so we will just describe the most basic processes involved in 
H2 formation. However, we note that the papers by Abel et al. (1997)^^ and 
by GP98 provide a much more accurate description of primordial chemistry. 

1.3.3.1 Atomic Hydrogen and free electrons 

Apart from molecule formation (see below), the main reactions involving 
Hydrogen (and Helium, to which we can apply all the arguments below) 
are ionizations and recombinations. 

^ The coUisional rate coefficients given in the Abel et al. (1997) paper can 
be readily obtained through a FORTRAN code available on the web page 
http://www.astro.psu.edu/users/tabel/PGas/LCA-CM.html ; also note that on Tom 
Abel's web site Ihttp:/ /www. tomabel.com/ I it is possible to find several useful informa- 
tions about the primordial universe in general. 
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Ionizations can be produced both by radiation + hv ^ H+ + e~) 
and by collisions, mainly with free electrons (H + e~ —> H"*" + 2e~) but 
also with other H atoms (H + H — > H+e~ +H). Photoionizations dominate 
over collisions as long as UV photons with energy above the H ionization 
threshold arc present (sec e.g. Ostcrbrock 1989), but this is not always the 
case before the formation of the first luminous objects, when only the CMB 
radiation is present; even after some sources of radiation have appeared, it 
is likely that their influence will be mainly local, at least until reionization. 
However, in the primordial universe the electron temperature is low, and 
coUisional ionizations are relatively rare. 

Recombinations (H+ +e~ — > H + hu) have much higher specific rates, 
since they do not have a high energy threshold. They probably dominate 



Table 1.1. Reaction rates for some of the most important reactions in primordial 
gas, plus the main reactions involved in the formation of HD. In the formulae, 
Tj is the temperature of the radiation field (for our purposes, the temperature 
of the CMB radiation) and temperatures need to be expressed in Kelvin; j(i^w) 
is the radiation flux (in ergs^^ cm^^) at the central frequency vtjw of the Ly- 
man- Werner bands, hi>Lw = 12.87 eV. The rates come from compilations given 
by Tegmark et al. 1997 (reactions 1-7), Palla et al. 1983 (reactions 8-13), Abel et 
al. 1997 (reaction 14) and Bromm et al. 2002 (reactions 15-19). 



Reaction 






Rate 




H+ -^e- 






fci ~ 


1.88 X io-^"T-""*'' cm^s"' 


H+e- 




R- + hu 


fca ^ 


1.83 X io~i«tO-88 cm^s-i 


H- +H 




H2 + e- 


fca ^ 


1.3 X 10"'^ cm'^s"^ 


H" +hu 




H + e" 


ki ~ 


0.114T2-i»e-«'"»/^-. cm-^s-i 


H+ +H 




Rt + hv 


fcs ^ 


1.85 X io-23t1« cm^s"^ 


H++H 




H2 - +H+ 


fce ^ 


6.4 X 10"^° cm^s"^ 


+ hi' 




H+ +H 


k, ~ 


6.36 X io5e-^i6oo/T, ^^^^3^-1 


H + H + H 




H2 + H 




5.5 X 10-29^-^ cmS"^ 


H2-hH 




H-hH-hH 


kg ~ 


6.5xlO-^r-i/2e-«2°°o/^x 










x(l-e-«°°°/^)cm8s-i 


H + H + H2 




H2 + H2 


kio =: 




H2 -|- H2 




H + H H2 


fell 


ife9/8 


H-l-e- 




U+ +e- + e- 


fel2 =^ 


i 5.8 X 10-"Ti/2e-is80oo/r cm^s-i 


H-l-H 




R+ + e- +R 


fel3 =^ 


i 1.7 X 10-*fci2 


H2 + hu 




H+ +H 


fcl4 =: 


i 1.1 X 10*j(fLw) S-1 


D+ + e- 




D + /w/ 


fel5 = 


; 8.4 X 10""T"°-^x 








x{ 




D + H+ 




D+ +H 


fcl6 =: 


i 3.7 X io-"'T»-2«e-«/^ cm^s-i 


D+ +H 




D + H+ 


fcl7 = 


i 3.7 X io-"'T"-2« cm»s-i 


D+ +H2 




H+ + HD 


fcl8 = 


; 2.1 X 10"'-' cm'^s"^ 


HD + H+ 




H2 + D 


fel9 = 


i 1.0 X lO-^e-W cm^s-i 
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the evolution of free electrons, which can be relatively abundant as residu- 
als from the recombination era {e.g. Peebles 1993), and are important for 
H2 formation (see below). 

Simple approximations for the reaction rates of coUisional ionizations 
and recombinations are given in Table [TTTl 



1.3.3.2 H2 formation and disruption 

At present, H2 is commonly believed to form mainly through reactions tak- 
ing place on the surface of dust grains (but see e.g. Cazaux & Spaans 2004). 
Such a mechanism cannot work in the primordial universe, when the metal 
(and dust) abundance was negligible. The first two mechanisms leading to 
formation of H2 in a primordial environment was described by McDowell 
(1961) and by Saslaw & Zipoy (1967); soon after the publication of this 
latter paper, H2 started to be included in theories of structure formation 
(such as the PD68 paper about globular cluster formation through H2 cool- 
ing). Both these mechanisms consist of two stages, involving either e~ or 
H+ as catalyzers. The first (and usually most important) goes through the 
reactions 

R + e-->ll-+hv (1.30) 
H-+H^H2-+e- (1.31) 

while the second proceeds as 

H++H^H+ + ;iiy (1.32) 
H++H^H2-+H+. (1.33) 

In both cases, H2 production can fail at the intermediate stage if a pho- 
todissociation occurs {R- + hiy ^ R + e-, or + hiy ^ R+ + H). The 
rates of all these reactions are listed in Table [TTTl 

By combining the reaction rates of these reactions, it is possible (see 
Tegmark et al. 1997; hereafter T97) to obtain an approximate evolution 
for the ionization fraction x = /n and the H2 fraction = 2nH2 /n 
(here n is the total density of protons, n ~ nn + + 2nH2): 

x{t) - (1.34) 

^ ' 1 + XQTlklt ^ ' 

/h,(0 -/o + 2^1n(l + a;onfcit) with (1.35) 
fci 

" k3 + ki/[n{l-x)] ke + kr/[n{l-x)] ^' ' 

where the various ki are the reaction rates given in Table [TTTl and xq and 
/o are the initial fractions of ionized atoms and of H2 molecules. 
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Another mechanism for H2 formation, which becomes important at 
(cosmofogicaUy) very high densities uh ^ 10^ cm" 3 are the so-cahed 3- 
body reactions described by Palla et al. (1983). While the previous mech- 
anisms are limited by the small abundance of free electrons (or of H"*"), 
reactions such as 

H + H + H^H2 + H (1.37) 

can rapidly convert all the hydrogen into molecular form, provided the 
density is high enough. For this reason, they are likely play an important 
role during the formation of a primordial protostar. 

Finally, H2 can be dissociated through coUisional processes such as 
reactions 9 in table 11.11 but probably the most important process is its 
photo-dissociation by photons in the Lyman- Werner bands (11.26-13.6 eV; 
but Abel et al. 1997 found that the most important range is between 12.24 
and 13.51 eV). These photons are below the H ionization threshold, there- 
fore they can diffuse to large distances from their source. So, any primordial 
source emitting a relevant number of UV photons {e.g. a ~ IOOMq star) 
is likely to have a major feedback effect, since it can strongly reduce the 
amount of H2 in the halo where it formed (and also in neighbouring halos), 
inhibiting further star formation. 

I. 3.3.3 Approximate predictions 

The above information can be used for making approximate predictions 
about the properties of the halos hosting the first luminous objects. Such 
kind of predictions was started by Couchman & Rees (1986), and especially 
by T97 and their example was followed and improved by several authors 
(see below). 

The basic idea is to have a simplified model for the evolution of the 
H2 fraction and the cooling rate inside spherical top- hat fluctuations, in 
order to check whether after virialization the baryons are able to cool (and 
keep collapsing), or they just "settle down" at about the virial radius. Such 
approximate models are fast to compute, and it is easy to explore the gas 
behaviour on a large range of virialization redshifts and fluctuation masses 
(or, equivalcntly, virial temperatures). In this way, for each virialization 
redshift it is possible to obtain the minimum molecular fraction which is 
needed for the collapse to proceed, and the minimum size of halos where 
such abundance is achieved. The results interestingly point out that the 
molecular fraction threshold separating collapsing and non-collapsing ob- 
jects has an almost redshift-independent value of fu^ 5 x lO"** (see fig. 

II. 12|l . Instead, the minimum halo mass actually evolves with redshift (see 

fig, rm . 

Predictions about the ability of the baryons inside each kind of halo 
to keep collapsing after virialization can then be combined with Press- 
Schechter predictions about the actual abundances of halos. For instance. 
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T [K] 

Figure 1.12. Comparison of tlie H2 fraction needed for an lialo to collapse 
and H2 fraction which can be formed inside the halo in an Hubble time (from 
T97). This is shown as a function of halo virial temperature and for three dif- 
ferent virialization redshifts {z = 100: solid; z = 50: short dashes; z = 25:long 
dashes). The three dots mark the minimum H2 abundance which is needed for 
collapse at the three considered redshift, and it can be seen that they all are at 
~ 5 X 10--*. 

in fig. 11.131 the solid black (almost vertical) line shows where the masses 
of 3cr fluctuations lie, as a function of redshift. So, if we decide to neglect 
the rare fluctuations at more than 3a from the average, that figure tells us 
that the first luminous objects can start forming only at z ;$ 30, in objects 
with a total mass ^ 2 x 10^ Mq. 

Such result is subject to a number of uncertainties, both about the 
"details" of the model and about the processes it neglects; here are some 
of the more interesting developments: 

(i) Abel et al. (1998) found that the minimum mass is strongly affected by 
the uncertainties in the adopted H2 cooling function, with differences 
that could reach a factor ~ 10 in the minimum mass estimate 

(ii) Fuller & Couchman (2000) used numerical simulations of single, high-cr 
density peaks in order to improve the spherical collapse approximation 
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Figure 1.13. Evolution with virialization redshift of the minimum mass which is 
able to cool, collapse and possibly form stars. Only the halos whose {zvIt, M) fall 
outside the filled region are able to cool fast enough. The filled region is made of 
a red part, where CMB radiation prevents the cooling, and a yellow part where 
the cooling is slow because of a dearth of H2 molecules. The two parallel dashed 
lines correspond to virial temperatures of 10'' K (the highest one) and 10^ K (the 
lowest one), while the almost vertical line in the middle of the plot corresponds 
to 3(7 peaks in a SCDM (S7m = 1, Q.a ~ 0) cosmology. 



(iii) Machacek et al. (2001) and Kitayama et al. (2001) investigated the 
influence of background radiation 

(iv) Yoshida et al. (2003) used larger-scale numerical simulations and found 
that also the merging history of an halo could play a role, since frequent 
mergings heat the gas and prevent or delay the collapse of ~ 30% of 
the halos. 

As a result of the improved modeling (and of a different set of cosmo- 
logical parameters as ACDM has substituted SCDM), in the most recent 
papers the value of the minimum halo mass for the formation of the first 
luminous objects is somewhat reduced to the range 0.5-1 xlO^ M©, with 
a weak dependence on redshift and a stronger dependence on other param- 
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eters (background radiation, merging history etc.). 

1.4 Numerical cosmological hydrodynamics 

Numerical simulations are an important tool for a large range of astro- 
physical problems. This is especially true for the study of primordial stars, 
given the absence of direct observational data about these stars, and the 
relatively scarce indirect evidence we have. 

1.4.1 Adaptive refinement codes (ENZO) 

The two problems which immediately emerge when setting up a simulation 
of primordial star formation are the dynamical range and the required 
accuracy in solving the hydrodynamical equations. 

When studying the formation of objects in a cosmological context, we 
need both to simulate a large enough volume (with a box size of at least 
100 comoving kpc), and to resolve objects of the size of a star (~ 10^^ cm 
in the case of the Sun), about 11 orders of magnitude smaller. 

This huge difference in the relevant scales of the problem is obviously 
a problem. It can be attenuated by the use of Smoothed Particle Hydro- 
dynamics (SPH; see e.g. Monaghan 1992), whose Lagrangian nature has 
some benign effects, as the simulated particles are likely to concentrate 
(i.e. , provide resolution) in the regions where the maximum resolution is 
needed. However, even if this kind of method can actually be employed (see 
Bromm et al. 1999, 2002), it has at least two important drawbacks. First 
of all, the positive effects mentioned above cannot bridge in a completely 
satisfactory way the extreme dynamical range we just mentioned, since the 
mass resolution is normally fixed once and for all at the beginning of the 
simulation. Second, SPH is known to have poor shock resolution proper- 
ties, which casts doubts on the results when the hydrodynamics becomes 
important. 

The best presently available solution for satisfying both requirements 
is the combination of an Eulerian approach (which is good for the hydrody- 
namical part) and an adaptive technique (which can extend the dynamical 
range). This is known as Adaptive Mesh Refinement (AMR; see e.g. Berger 
& Colella 1989; Norman 2004), and basically consists in having the simula- 
tion volume represented by a hierarchy of nested grids (meshes) which are 
created according to resolution needs. 

In particular, the simulations we are going to describe in the following 
paragraphs were made with the code ENZO^ (see O'Shea et al. 2004 and 
references therein for a full description) . Briefly, ENZO includes the treat- 
ment of gravitational physics through N-body techniques, and the treat- 

^ The ENZO code can be retrieved at the web site |http: / / cosmos .ucsd.edu/ enzo / 1 
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Figure 1.14. Overview of tlie evolution leading to the formation of a primordial 
star. The top row shows the gas density, centered at the pre-galactic object within 
which the star is formed. The four projections are labelled with their redshifts. 
Pre-galactic objects form from very small density fluctuations and continuously 
merge to form larger objects. The middle and bottom rows show thin slices 
through the gas density and temperature at the final simulation stage. The four 
pairs of slices are labelled with the scale on which they were taken, starting from 
6 (proper) kpc (the size of the simulated volume) and zooming in down to 0.06 pc 
(12,000 AU). In the left panels, the larger scale structures of filaments and sheets 
are seen. At their intersections, a pre-galactic object of ~ 10^ Mq is formed. The 
temperature slice (second panel, bottom row) shows how the gas shock heats as it 
falls into the pre-galactic object. After passing the accretion shock, the material 
forms hydrogen molecules and starts to cool. The cooling material accumulates 
at the centre of the object and forms the high-redshift analog to a molecular 
cloud (third panel from the right), which is dense and cold (T ~ 200K). Deep 
within the molecular cloud, a core of ~ 100 Mq, a few hundred K warmer, is 
formed (right panel) within which a ~ 1 Mq fully molecular object is formed 
(yellow region in the right panel of the middle row), (from ABN02). 
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ment of hydrodynamics through the piecewise parabohc method (PPM) 
of Woodward & ColeUa 1984 (as modified by Bryan et al. 1995 in order 
to adapt to cosmological simulations). ENZO can optionally include the 
treatment of gas cooling (in particular H2 cooling, as described in the 
preceeding sections), of primordial non-equilibrium chemistry (see the pre- 
ceeding sections) and of UV background models {e.g. the ones by Haardt 
& Madau 1996) and heuristic prescriptions (see Cen & Ostriker 1992) for 
star formation in larger-scale simulations. 

1.4.2 Formation of the first star 

The use of AMR codes for cosmological simulations of the formation of 
the first objects in the universe was pioneered by Abel et al. (2000) and 
further refined by Abel et al. (2002; hereafter ABN02), where the dy- 
namic range covered by the simulations was larger than 5 orders of mag- 
nitude in scale length, i.e. the wide range between an almost cosmological 
( ^ 100 comoving kpc, i.e. ^5 proper kpc) and an almost stellar scale 
( ^ 1000 AU ~ 10"^ pc). In the whole process, the AMR code kept intro- 
ducing new (finer resolution) meshes whenever the density exceeded some 
thresholds, or the Jeans length was resolved by less than 64 grid cells. 

The simulations were started at a redshift z = 100 from cosmolog- 
ically consistent initial conditions®, and the code also followed the non- 
equilibrium chemistry of the primordial gas, and included an optically thin 
treatment of radiative losses from atomic and molecular lines, and from 
Compton cooling. 

The main limitation of these simulations was the assumption that the 
cooling proceeds in the optically thin limit; such assumption breaks down 
when the optical depth inside H2 lines reaches unity (corresponding to a 
Jeans length ^ 10'^ AU ^ 0.01 pc). However, the simulations were halted 
only when the optical depth at line centres becomes larger than 10 (Jeans 
length of about - 10 AU (~ 10"'' pc), since it was unclear whether Doppler 
shifts could delay the transition to the optically thick regime. 

1.4.2.1 Summary of the evolution: radial profiles 

In Figures 11.141 and 11.151 we show the evolution of gas properties both in 
pictures and in plots of spherically averaged quantities, as presented in Abel 
et al. (2002). From these figures, and in particular from the local minima 
in the infall velocity (Fig. 11.15b ) it is possible to identify four characteristic 
mass scales: 

^ Such conditions were taken from an SCDM model with fi^ = 0, Qm = li f^b = 0.06, 
-ffo = 50 km/sMpc~^ which is quite different from the "concordance" ACDM model. 
However, the final results are believed to be only marginally affected by differences in 
these cosmological parameters. 
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Figure 1.15. Radial mass-weighted averages of physical quantities at seven 
different simulation times. (A) Particle number density in cm~^ as a function of 
radius; the bottom line corresponds to 2 = 19, and moving upwards the "steps" 
from one line to the next are of 9 x 10** yr, 3 x 10^ yr, 3 x lO" yr, 3000 yr, 1500 yr, 
and 200 yr, respectively; the uppermost line shows the simulation final state, at 
z = 18.181164. The two lines between 0.01 and 200 pc give the DM mass density 
(in GeVcm""^) at 2; = 19 and the final time, respectively. (B) Enclosed gas mass. 
(C) Mass fractions of H and H2 . (D) Temperature. (E) Radial velocity of the 
baryons; the bottom line in (E) shows the negative value of the local speed of 
sound at the final time. In all panels the same output times correspond to the 
same line styles, (from ABN02). 
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(i) The mass scale of the pre-galactic halo, of 7 x 10^ Mq in total, 
consistent with the approximate predictions discussed in the previous 
sections 

(ii) The mass scale of a "primordial molecular cloud", ~ 4000Mo: the 
molecular fraction in this region is actually very low ( ^ 10"'^), but it is 
enough to reduce the gas temperature from the virial value ( 10"^ K) 
to - 200 K 

(iii) The mass scale of a "fragment", ^ IOOA/q, which is determined by 
the change in the H2 cooling properties at a density n ^ cm~^ 
(for a complete discussion, see Bromm et al. 1999, 2002), when Local 
Thermal Equilibrium is reached and the cooling rate dependence on 
density flattens to A (x n (from A oc n^); this is also the first mass 
scale where the gas mass exceeds the Bonnor-Ebert mass (Ebert 1955, 
Bonnor 1956) MBE{T,n,fi) ~ eiMoT^/^n-i/^^-^ (with T in Kelvin 
and n in cm~^), indicating an unstable collapse 

(iv) The mass scale of the "molecular core", ~ IM© is determined by 
the onset of 3-body reactions at densities in the range n ^ 10^ — 
10^" cm"'', which leads to the complete conversion of hydrogen into 
molecular form; at this stage, the infall flow becomes supersonic (which 
is the precondition for the appearance of an accretion shock and of a 
central hydrostatic flow); the increase in the H2 abundance due to the 
formation of this molecular core also leads to the transition to the 
optically thick regime. 

1.4-2.2 Angular momentum 

In Fig. 11.161 we show the evolution of the radial distribution of average 
specific angular momentum (and related quantities). It is remarkable that 
in ABN02 simulations rotational support does not halt the collapse at any 
stage, even if this could be a natural expectation. 

There are two reasons for such (apparently odd) fact: 

(i) As can be seen in panel A of Fig. 11.161 the collapse starts with the 
central gas having much less specific angular momentum than the av- 
erage (this is typical of halos produced by gravitational collapse; see 
e.g. Quinn & Zurek 1988). That is, the gas in the central regions starts 
with little angular momentum to lose in the first place. 

(ii) Second, some form of angular momentum transport is clearly active, 
as is demonstrated by the decrease in the central specific angular mo- 
mentum. Turbulence is likely to be the explanation: at any radius, 
there will be both low and high angular momentum material, and re- 
distribution will happen because of pressure forces or shock waves: 
lower angular momentum material will selectively sink inwards, dis- 
placing higher angular momentum gas. It is notable that this kind of 
transport will be suppressed in situations in which the collapse occurs 




Figure 1.16. Radial mass weigiited averages of angular momentum-related 
quantities at different times (the same as in fig. If .f 5^ . (A) specific an- 
gular momentum L. (B) Rotational speed in units of Keplerian velocity 
VKc-p = (GMr/r)^/^. (C) Rotational speed [L/r). (D) Rotational speed in units 
of the sound speed Cs (from ABN02). 



on the dynamical time scale, rather than the longer cooling time scale: 
this is the likely reason why this kind of mechanism has not been ob- 
served in simulations of present day star formation (see e.g. . Burkert 
& Bodenheimer 1993). 



1.4.3 SPH results 



Bromm et al. (1999, 2002) have performed simulations of primordial star 
formation using an SPH code which included essentially the same physical 
processes as the simulations we just described. Apart from the numerical 
technique, the main differences were that they included deuterium in some 
of their models, and that their initial conditions were not fully cosmological 
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{e.g. , since they choose to study isolated halos, the angular momenta are 
assigned as initial conditions, rather than generated by tidal interactions 
with neighbouring halo). 

In Fig. I1.17l and ll.l8l we show the results of one of their simulations, 
which give results in essential agreement with those we have previously 
discussed. 

An interesting extra-result of this kind of simulation is that the authors 
are able to assess the mass evolution of the various gaseous clumps in the 
hypothesis that feedback is unimportant (which could be the case if the 
clumps directly form intermediate mass black holes without emitting much 
radiation, or if fragmentation to a quite small stellar mass scale happens). 
They find that a significant fraction (~ 0.5) of the halo gas should end 
up inside one of the gas clumps, although it is not clear at all whether 
this gas will form stars (or some other kind of object). Furthermore, they 
find that clumps are likely to increase their mass on a timescale of about 
lO"^ yr (roughly corresponding to the initial time scale of the simulated 
halo), both because of gas accretion and of mergers, and they could easily 
reach masses <L 10^ Mq. Obviously, this result is heavily dependent on 
the not very realistic assumed lack of feedback. 

1.5 Protostar formation and accretion 

Full three-dimensional simulations (such as the ones of ABN02) are not 
able to reach the stage when a star is really formed. They are usually 
stopped at densities (n ~ 10^" — 10^^ cm~'^) which are much lower than 
typical stellar densities {p ^ 1 gcm"'^, n ~ 10^^ cm^^). In fact, at low 
densities the gas is optically thin at all frequencies, and it is not necessary 
to include radiative transfer in order to estimate the gas cooling. Instead, 
at densities n ^ 10^*^ cm~'^ some of the dominant H2 roto-vibrational lines 
become optically thick and require the treatment of radiative transfer. 

For this kind of problem, this is a prohibitive computational burden, 
and at present the actual formation of a protostar can not be fully investi- 
gated through self-consistent 3-D simulations. In order to proceed further, 
it is necessary to introduce some kind of simplification in the problem. 

1.5.1 Analytical results 

Historically, there were several studies based on analytical arguments {i.e. 
stability analysis) or single-zone models, leading to different conclusions 
about the properties of the final object. 

An early example is the paper by Yoneyama (1972), in which it was ar- 
gued that fragmentation takes place until the opacity limit is reached. How- 
ever, Yoneyama looked at the opacity limit of the entire "cloud" (roughly 
corresponding to one of the 10^ — 10^ Mq mini-halos we consider at present; 
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Figure 1.17. Typical result of the SPH simulations by Bromm et al. 2002. The 
figure shows the morphology of a simulated 2 x 10® M© halo virializing at a ~ 30 
just after the formation of the first clump of mass 1400 Mq (which is likely 
to produce stars). The two top row panels shows the distribution of the dark 
matter (which is undergoing violent relaxation). The two bottom panels show 
the distribution of gas, which has developed a lumpy morphology and settled at 
the centre of the potential well. 



originally, this was mass the scale proposed by PD68) rather than the pu- 
tative fragments and arrived at the conclusion that fragmentation stopped 
for masses ^ QOMq. 

More recently, Palla et al. 1983 pointed to the increase of the cooling 
rate at n ^ 10* cm~^ (due to H2 fast formation through 3-body reactions) 
as a possible trigger for instability, leading to fragmentation on mass scales 



Protostar formation and accretion 35 



10 




-2 02458 -2 02468 

log n„ [crm"^] log n„ [cm"^] 



Figure 1.18. Gas properties in the SPH simulations by Bromm et al. 2002. The 
figure shows the properties of the particles shown in Fig. 11.171 Panel (a) shows 
the free electron abundance X Q £LS 3, fraction of H number density. Panel (b) 
shows the H2 abundance /hj (note that even the highest density particles have 
nn ~ 10* cm^, so 3-body reactions are unimportant and /hj ~ 10~^). Panel (c) 
shows the gas temperature (low density gas gets to high temperatures because 
H2 cooling is inefficient). Panel (d) shows the value of the Jeans mass which can 
be obtained by using each particle temperature and density. All the panels have 
the hydrogen number density UrmH on the X axis. 



of 0.1 M©; such instability has actually been observed in the simulations 
of Abel et al. (2003), but it does not lead to fragmentation because its 
growth is too slow (see also Sabano & Yoshii 1977, Silk 1983, Omukai 
& Yoshii 2003 and Ripamonti & Abel 2004 for analytical fragmentation 
criteria and their application to this case). 
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Figure 1.19. Proto-stellar collapse, formation of an hydrostatic core and start 
of the proto-stellar accretion phase as can be found with 1-D simulations. The 
four panel show the profiles of density (top left; a), infall velocity (top right; 
b), temperature (bottom left;d) and H2 abundance (bottom right; d) as a func- 
tion of enclosed mass at 10 different evolutionary stages (0=initial conditions; 
9=final stage of the computation). The most relevant phases are the rapid for- 
mation of H2 (2-3) and the formation of a shock on the surface of the hydrostatic 
core (5-6-7), followed by the onset by accretion. Also note the almost perfect 
power-law behaviour of the density profile before core formation, (from Ripar 
monti et al. 2002). 




Figure 1.20. Comparison of tiie accretion and tiie Kelvin-Helmoltz time 
scales for a primordial protostars, as a function of protostellar mass. The 
Kelvin-Helmoltz contraction time (obtained using the ZAMS luminosity as given 
by Schaerer 2002) is shown as the solid black line, while the dashed line and the 
solid line with circles show the time which is needed to accrete each mass of gas; 
they are based on the results of ABN02 and differ slightly in the way in which 
they were obtained. The dotted lines mark three constant accretion rates of 10~^ 
(bottom line), 5 x 10~^ (middle line) and 10~^ Mo/yr (top line). 



1.5.2 Mono-dimensional models 

1.5.2.1 The formation of the protostellar core 

A less simplified approach relies on ID studies, such as those of Omukai & 
Nishi (1998) and Ripamonti et al. (2002): in such studies the heavy price of 
assuming that the collapsing object always remains spherical (which makes 
impossible to investigate the effects of angular momentum, and prevents 
fragmentation) is compensated by the ability to properly treat all the other 
physical processes which are believed to be important: first of all, the de- 
tailed radiative transfer of radiation emitted in the numerous H2 resonance 
lines, then the transition to continuum cooling (at first from molecular 
Collision-Induced Emission, and later from atomic processes) and finally 
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Figure 1.21. Evolution of the protostellar radius as a function of protostellar 
mass in the models of Omukai & Palla (2003). The models differ in the assumed 
accretion rates. The dotted curves show the results of models where A/corc is 
assumed to be constant at values of 1.1 x 10~^ M0/yr (dotted curve starting at 
i?corc ^ 30 Rq), 2.2x10"^ Mo/yr (starting at i?corc ^ 40 Rq), 4.4 x 10"^ Mg/yr 
("fiducial" model; starting at -Rcore — 50 Rq) and 8.8 x 10"'^ Mg/yr (starting at 
_Rcorc — 65 Rq). The solid thick curve shows the evolution of a model in which 
the accretion rate was obtained from the extrapolation of the ABN02 data. 



the non-ideal equation of state which becomes necessary at high densities, 
when the gas becomes increasingly hot and largely ionized. 

Such studies find that the collapse initially proceeds in a self-similar 
fashion (in good agreement with the solution found by Larson 1969 and 
Penston 1969); at the start of this phase (stages 1 and 2 in fig. 11.191 to 
be compared with the profiles at corresponding densities in fig. I1.15|) their 
results can be compared with those of 3-D simulations, and the agreement 
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is good despite the difference in the assumed initial conditions: this is likely 
due to the self-similar properties. At later stages the comparison is not as 
good, but this is presumably due to the differences in the cooling rate (the 
3D models at n ^ 10^° cm^'^ only include the optically thin cooling rate). 

It is noticeable that the transition to optically thick cooling does not 
stop the self-similar phase, contrary to theoretical expectations (see e.g. 
Rees 1976). The reason is that the reduction in the cooling rate is smooth, 
and is further delayed by the onset of CIE cooling and by the start of H2 dis- 
sociation, which acts as an heat well (each H2 dissociation absorbs 4.48 eV 
of thermal energy) and provides an alternative "disposal" mechanism for 
the excess thermal energy. 

The self-similar phase proceeds until the central density is n ^ 10^^ cm~ 
^ 10^^ gcm^'^, when a small (Mcoro,o ~ 3 x 10~^ Mq) hydrostatic core 
eventually forms. Such core starts accreting the surrounding gas at a very 
high rate (Mcore ~ 0.01 — 0.1 AfQyr"^), comparable to the theoretical 
predictions of Stabler et al. 1980 

M_ - cl/G . 4 X 10-3 (^) (^) yr- (1.38) 

where Cs = [k-QT/ (^m-a)]^^'^ is the isothermal sound speed of a gas with 
mean molecular weight fi. 

Although the Omukai & Nishi (1998) and the Ripamonti et al. (2002) 
studies predict that such a high accretion rate will decline (approximately 
as Mcoro oc t^'^-^), it is clear that, if the accretion proceeds unimpeded, it 
could lead to the formation of a star with a mass comparable to the whole 
proto-stellar cloud where the proto-star is born, that is, about lOOM©; such 
a process would take a quite short time of about 10^ — 10^ yr. 

This scenario could be deeply affected by feedback effects: even not 
considering the radiation coming from the interior of the proto-star, the 
energy released by the shocked material accreting onto the surface can 
reach very high values, comparable with the Eddington luminosity, and is 
likely to affect the accretion flow. 

Figure [1.201 compares the accretion time scale (which can be extrap- 
olated from the data at the end of the ABN02 simulations) to the Kelvin- 
Helmoltz timescale, 

tKK = ^ (1.39) 

where L is the luminosity of the protostar. This plot tells us that the 
accretion timescale is so fast that the stellar interior has very little space for 
a "re-adjustment" (which could possibly stop the accretion) before reaching 
a mass of - 10 - lOOMg. 

However, it can be argued that such a readjustment is not necessary, 
since in the first stages most of the protostellar luminosity comes from the 
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accretion process itself: 



where we have inserted reahstic values for Mcore, Mcorc and i?corc- This 
luminosity is quite close to the Eddington luminosity (LEdd — 1-3 x 10'^* 
(M/Mq) ergs~^), and radiation pressure could have a major effect. 

Ripamonti et al. (2002) found that this luminosity is not able to stop 
the accretion, but they do not properly trace the internal structure of the 
core, so their results cannot be trusted except in the very initial stages of 
accretion (say, when Mcore ~ O.IM©). 

A better suited approach was used in studies by Stahler et al. (1986) 
and, more recently, by Omukai & Palla (2001, 2003), who assumed that 
the accretion can be described as a sequence of steady-state accretion flows 
onto a growing core. The core is modelled hydrostatically, as a normal 
stellar interior (including the treatment of nuclear burning), while the ac- 
creting envelope is modelled through the steady state assumption, in con- 
junction with the condition that outside the "photosphere" (the region 
where the optical depth for a photon to escape the cloud is ^1) the ac- 
creting material is in free-fall. As shown in Fig. 11.211 Omukai & Palla 
(2003) find that even if feedback luminosity deeply affects the structure 
of the accreting envelope, it is never able to stop the accretion before the 
proto-stellar mass reaches 60 — 100 Mq as a minimum; after that, the final 
mass of the protostar depends on the assumed mass accretion rate Mcore ^ 
for Mcore ^ 4 X 10^^ Mq yr~^, the accretion can proceed unimpeded until 
the gas is completely depleted (or the star explodes as a supernova); oth- 
erwise, the radiation pressure is finally able to revert the accretion; with 
high accretion rates this happens sooner, and the final stellar mass is rela- 
tively low, while for accretion rates only slightly above the critical value of 
~ 4 X 10""^ Mq yr~^ the stellar mass can reach about 300 Mq. 

If such predictions are correct, the primordial Initial Mass Function 
is likely to be much different from the present one, reflecting mainly the 
mass spectrum of the gas fragments from which the stars originate, and a 
mass of ~ 100 Mq could be typical for a primordial star. However, we note 
that this important results could change as a result of better modeling. For 
example, deep modifications of the envelope structure are quite possible if 
a frequency-dependent opacity (rather than the mean opacity used in the 
cited studies) is included. 
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1.6 Discussion 

The previous sections show that, although not certain at aU (because of 
the big uncertainty about feedback effects), numerical simulations tend to 
favour the hypothesis that primordial stars had a larger typical mass than 
present-day stars. 

If this is true, it could indeed solve some observational puzzle, such 
as why wc have never observed a single zcro-mctallicity star (answer; they 
have exploded as supernovae and/or transformed into compact objects at 
a very high redshift), and maybe help explaining the relatively high metal- 
licities {Z ^ 10~^Zq even in low column density systems) measured in the 
Lyman a forest (see e.g. Cowie & Songaila 1998), or the high ('^ solar) 
metallicities we observe in the spectra of some quasars already at redshift 
6 (Fan et al. 2000, 2001, 2003, 2004; Walter et al. 2003). 

However, a top-heavy primordial IMF also runs into a series of prob- 
lems, which we will discuss in the remaining of this section. 

1.6.1 UV Radiation feedback 

First of all, if a moderately massive star (even M ^ 20— 3OM0 is likely to be 
enough) actually forms in the centre of an halo in the way shown by ABN02, 
it will emit copious amounts of UV radiation, which will produce important 
feedback effects. Indeed, the scarcity of metals in stellar atmospheres (see 
Schaerer 2002, and also Tumlinson & ShuU 2000, Bromm Kudritzki & Loeb 
2001) results in UV fluxes which are significantly larger than for stars of the 
same mass but with "normal" metallicity. Furthermore, this same scarcity 
of metals is likely to result in a negligibly small density of dust particles, 
further advancing the UV fluxes at large distances from the sources when 
compared to the present-day situation. 

Massive primordial objects are also likely to end up into black holes 
(either directly or after having formed a star^), which could be the seeds of 
present day super-massive black holes (see 6.17. Volonteri, Haardt & Madau 
2003). If accretion can take place onto these black holes (also known as 
mini- quasars), they are likely to emit an important amount of radiation 
in the UV and the X-rays (this could also be important in explaining the 
WMAP result about Thomson scattering optical depth; see e.g. Madau et 
al. 2003). 

UV photons c;an have a series of different effects. First of all, we al- 
ready mentioned in the section about chemistry that Lyman- Werner pho- 
tons (11.26 eV < hv < 13.6 eV) can dissociate H2 molecules, preventing 

^ Here we think of a star as an object where quasi-hydrostatic equilibrium is reached 
because of stable nuclear burning; according to this definition, a blaclc liolc can be formed 
without passing through a truly stellar phase, even if it is very likely to emit significant 
amounts of radiation in the process. 
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Figure 1.22. Average flux (in units of ergcm^^ sr^^ Hz^^) at an observation 
redsliift 2;obs = 25 coming from sources turning on at z = 35. Tlie top panel shows 
the effects of absorption of neutral hydrogen and helium, strongly reducing the 
flux between ~ 13.6 eV and ~ 1 keV (the solid and dashed lines show the 
absorbed and unabsorbed flux, respectively). The bottom panel shows the same 
quantities in a much smaller energy range, in order to illustrate the sawtooth 
modulation due to line absorption below 13.6 eV. (from Haiman, Rees & Loeb 
1997). 



further star formation. Such photons are effective even at large distances 
from their sources, and even in a neutral universe, since their energy is 
below the H ionization threshold; the only obstacle they find on their way 
are some of the Lyman transitions of H, which can scatter these photons or 
remove them from the band; this results in a "sawtooth" modulation (see 
fig- ll.22|) of the spectrum. Haiman, Rees & Loeb (1997) argue that this 
negative feedback could conceivably bring primordial star formation to an 
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Figure 1.23. Dynamical evolution of the HII region produced by a 200 Mq star. 
Panels refer to ionization fraction profiles (top left), temperature distributions 
(top right), densities (bottom left) and velocities (bottom right). The dashed 
line in the density panel is for the Stromgren density (the density required to 
form a Stromgren sphere and therefore initially bind the ionization front) at a 
given radius. Profiles are given at times (from top to bottom in the density 
panel; from left to right in the others) (density panel), 63 kyr (all panels), 82 
kyr (ionization and temperature panels), 95 kyr (ionization panel) 126.9 kyr (all 
panels), 317 kyr (all panels), 1.1 Myr (temperature, density and velocity panels) 
and 2.2 Myr (all panels), which is the approximate main sequence lifetime of a 
200 Mq star. (From Whalen, Abel & Norman 2004). 

early stop; however, other authors found that this is not necessarily the 
case (Ciardi et al. 2000; Omukai 2001; Glover & Brand 2001), or that the 
feedback effect on H2 abundance can be moderated by the effects of X-ray 
photons coming from mini-quasars (Haiman Abel and Rees 2000; Glover 
& Brand 2003). 

A second obvious effect from UV photons is to ionize the interstellar 
and the intergalactic medium; it is unclear whether these kind of objects are 
important in the reionization history of the universe, but they will definitely 
form HII regions in their immediate neighbourhoods. In fig. II. 231 we show 



Discussion 44 



the results of Whalen, Abel & Norman 2004 about the evolution of an HII 
region produced by a 2OOM0 star inside an halo with the same density 
profile as found in the ABN02 paper: in the beginning the ionization front 
is D-type (that is, it is "trapped" because of the density) and expands 
because a shock forms in front of it; later (after about 70 kyr from the 
start of the UV emission) the ionization front reaches regions of lower 
densities and becomes R-type (radiation driven), expanding much faster 
than the shock. However, the shock keeps moving at speeds of the order 
of 30 kms~^ even after the source of ionizing photons is turned off. Such a 
speed is much larger than the rotational velocities of the mini-halos where 
primordial star formation is supposed to be taking place ( ^ 10 kms^^), so 
this shocks are likely to expel a very large fraction of the original gas content 
of the mini-halo. UV emission could even lead to the photo- evaporation 
of neighbouring mini-halos, similar to what Barkana & Loeb 1999 (see 
also Shapiro, Iliev & Raga 2004) found in the slightly different context to 
cosmological reionization. 

Since the dynamical timescale of a mini-halo is ^ 10^ yr and is longer 
than the timescale for this kind of phenomena (definitely shorter than the 
^1 — 10 Myr main sequence lifespan of massive stars, and very likely to 
be of the order of ~ 10"* — 10^ yr), the star formation in one mini- halo 
is likely to stop almost immediately after the first (massive) object forms. 
This means that, unless two or more stars form exactly at the same time 
(within ~ 1% of the mini-halo dynamical time), each mini- halo will form 
exactly one massive star. 

1.6.2 Supernovae feedback and metallicities 

After producing plenty of UV photons during their life, primordial massive 
stars are likely to explode as supernovae. This must be true for some 
fraction of them, otherwise it would be impossible to pollute the IGM with 
metals, as mass loss from zero-metal stars is believed to be unimportant 
(this applies to stars with M ^ 500 Mq; see Baraffe, Heger & Woosley 
2001, and also Heger et al. 2002); however, it is clearly possible that some 
fraction of primordial stars directly evolve into black holes. 

In figure 11.241 we show the results about the fate of zero metallicity 
stars, as obtained by Heger & Woosley 2002 (but see also Heger et al. 2003), 
indirectly confirming this picture and suggesting that pair-instability su- 
pernovae could play a major role if the primordial IMF is really skewed 
towards masses ^ 100 Mq. 

This has a series of consequences. First of all, supernova explosion are 
one of the very few phenomena directly involving primordial stars which we 
can realistically hope to observe. Wise & Abel (2004; see also Marri & Fer- 
rara 1998 and Marri, Ferrara & Pozzetti 2000 for the effects of gravitational 
lensing) investigated the expected number of such supernovae by means of 
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Figure 1.24. The "final mass function'' of non-rotating metal-free stars as a 
function of their initial mass, as given by from Heger & Woosley 2002. The 
thick black curve gives the final mass of the collapsed remnant, and the thick 
gray curve and the mass of the star at the start of the event that produces that 
remnant (mass loss, SN etc.); for zoro-mctal stars this happens to be mostly 
coincident with the dotted line, corresponding to no mass loss. Below initial 
masses of 5 — 10 AIq white dwarfs arc formed; above that, initial masses of up 
to ~ 25 Mq lead to neutron stars. At larger masses, black holes can form in 
different ways (through fall-back of SN ejecta or directly). Pair instability starts 
to appear above ~ 100 Mq, and becomes very important at ~ 140 Mq: stars 
with initial masses in the 140 — 260 Mq range are believed to completely disrupt 
the star, leaving no remnant; instead, above this range pair instability is believed 
to lead to a complete collapse of the star into a black hole. 
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a semi-analytic model combining Press-Schechter formalism, an evolving 
minimum mass for star forming halos and negative feedbacks, finding the 
rates shown in Fig. If .251 if such objects are pair-instability supernovae 
with masses ^ f 75 Mq, they should be detectable by future missions such 
as JWST*; if some of them are associated to Gamma Ray Bursts, the re- 
cently launched SwiflP satellite has a slim chance of observing them at a 
redshifts < 30 (see Gou et al. 2004). 

Supernova explosions obviously have hydrodynamical effects on the 
gas of the mini-halo where they presumably formed, especially in the case 
of the particularly violent pair-instability supernovae (see e.g. Ober, El Elid 
& Fricke 1983), and they are likely to expel it even if it was not removed by 
the effects of UV radiation (see e.g. MacLow & Ferrara 1999 and Ferrara 
& Tolstoy 2000), further reducing the probability of having more than 1 
star per mini-halo. The similarity with UV effects extends to the influence 
on neighbouring halos which can be wiped out by SN explosions (Sigward, 
Ferrara & Scannapieco 2004) . 

Finally, supernovae provide a mechanism for spreading metals into the 
IGM, as discussed e.g. by Madau, Ferrara & Rees 2001. In turn, these met- 
als will modify the cooling properties of the gas, and lead to star formation 
with a normal (Salpeter-like) IMF. If this is true, and if pair-instability 
supernovae dominate the yields from primordial stars, we should be able 
to distinguish the peculiar pair-supernova abundance pattern (described in 
Heger & Woosley 2001) as a "signature" of primordial origin. 

Searches for low-metallicity stars have a long history (see e.g. Beers 
1999, 2000), and their inability to find stars with mctallicities Z < 10~^Zq 
led to speculation that this metallicity marked the transition to a nor- 
mal IMF (see e.g. Bromm et al. 2001, Schneider et al. 2002). However, 
Christlieb et al. 2002 finally reported the discovery of an extremely metal 
poor star (HE0107-5240) with [Fe/H] = -5.3^°, a level compatible with 
a "second generation" star (i.e. , a star formed from material enriched by 
very few supernovae; for comparison. Wise & Abel 2004 find that primordial 
SNe could enrich the IGM to [Fe/H] ^ —4.1, although this is probably 
an upper limit). The abundance patterns in this star are quite strange 
(for example, the carbon abundance is slightly less than 1/10 solar), but a 
much better fit is provided by supernova yields of moderately massive stars 
(15 — 40 Mq) rather than from yields predicted for pair-instability super- 
novae coming from very massive stars. However, at the moment no model 
can satisfactorily fit all the observed abundances (see Bessel, Christlieb & 
Gustafsson 2004). 

Even if the primordial nature of this star must still be established, it 
represents a cautionary tale about the currently favoured predictions of a 

^ James Webb Space Telescope; more information at |http:/ /www.jwst.nasa.gov| 
^ More information at http://swift.gsfc.nasa.gov 
10 [Fe/H] = log.oiNpJNa) - log,o{NFe/NH)e- 
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Figure 1.25. Primordial supernova properties as reported by Wise & Abel 2004. 
The right panels show the diflterential rate of primordial supernovae per unit 
redshift (top) and the cumulative rate (bottom); both are per year of observation 
and per square degree. The left panels show the critical halo mass (in A/0) 
for primordial star formation (top), the comoving number density (in Mpc~^) 
of halos above the critical mass (middle) and the predicted specific intensity of 
the soft UV background in the Lyman-'Werner band (in erg s~^ cm~^ Hz~^ sr~^, 
bottom). The three lines in each panel refer to fixed primordial stellar masses of 
IOOM0 (solid), 2OOA/0 (dotted) and SOOMq (dashed). 

large number of massive or very massive primordial stars, and a reminder 
that better theoretical modeling is still needed, with particular regard to 
feedback effects during the stellar accretion phase. 
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